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Abstract 

The partial transposition(PT) operation is an effecient tool in detecting the 
inseparability of a mixed state. We give an explicit formula for the FT op- 
eration for the continuous variable states in Fock space. We then give the 
necessary and sufficient condition for the positivity of Gaussian operators. 
Based on this, a number of creterions on the inseparability and distillabil- 
ity for the multimode Gaussian states are naturally drawn. We finally give 
an explicit formula for the state in a subspace of a global Gaussian state. 
This formula, together with the known results for Gaussian states, gives the 
criterions for the inseparability and distillability in a subspace of the global 
Gaussian state. 

Quantum entanglement or inseparability plays a fundamentally important role in the fasci- 
nating area of quantum computation and quantum information. To check the inseparability 
of a (mixed) state, a particularly elegant mathematical tool, partial transposition(PT) op- 
eration was given by Peres |]l|. This method was later on applied to various discrete states 
0-^. While a lot of studies in quantum computation and quantum information(QCI) are 
based on the discrete states, recently a number of new proposals in the area were raised by 
using the continuous variable states In certain respects, the continuous variable states 

can have advantages to the discrete states. Consequently, the entanglement detection for 
continuous variable states is an interesting and important topic. The distillability and and 
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inseparability for Gaussian states were broadly studied |1T0|-|T^,P^ recently. In these works, 



the starting point is the Wigner function, whose PT operation is simple and clear. Actually, 
the PT operation in Fock space is also simple. In this letter, we directly use Peres' PT 
operation to the density operator itself. We shall first give an explicit formula for the PT 
operation for the continuous variable states in Fock space. We then check the inseparability 
of various types of practically existing states. We formulate a new type of practical impure 
squeezed state in the form of mixture of different squeezed vacuum states and Fock number 
states. We show this state and the squeezed state in phase damping channel are always 
inseparable. We then apply the PT operation to the Gaussian operator. We shall give the 
general criterion on the positivity of multi-variable Gaussian states in Fock space. This can 
be used to give the necessary and sufficient condition for inseparability or distillability for 
N ^ 1 system |T^ or arbitrary bipartite system []r3| , respectively. More generally, with the 



positivity condition for Gaussian state in Fock space we can easily investigate the entangle- 
ment prpoerty in a subspace of a global Gaussian state. 

The PT operation in Fock space. We start from the number state \ni) in Fock space, i 
from 1 to /. Without loss of any generality, operators in Fock space can be described by the 
summation of /-mode ket bra operator 

P=Y1 C{n„n'j\nun2 ■ ■ ■ni){n[,n'2 ■■■n'i\ + h.c. (1) 

{rii} 

Suppose mode 1 to j — 1 belong to subspace A and mode j to / belong to subspace B. The 
PT operation to subspace B is 

^PT ^ j2 C{n„n'j\ni,n2 " ■ ■ Uj^i, Uj, u'j^^ ■ ■ ■ u'l) {n[, u'^ ■ ■ ■n'j_^,nj,nj+i, ■ ■ ■,n;| + h.c. (2) 

{rii} 

An operator in Fock space can also be expressed in creation(a"'') and annihilation(a) oper- 
ators. Considering the fact = J2'i^=o + + and a = I]^o + l\n){n + 1|, 
we have a = a' ; a' = a , where k is an arbitrary natural number, the superscript T 
indicates the transpose operation. Without any loss of generality, an /-mode density oper- 
ator in Fock space can be written in the following normally ordered form (i.e., all creation 
operators in the left and all annihilation operators in the right) 

W = f{a\, ai, 4, 02 ■ ■ ■ aj, ai) = ^ c^i,,,,...^,,!.,^'''^! ■ ■ ■ al'''a'{' + h.c. (3) 

Note that the operators in the different mode are commute to each other. If we take the PT 
operation to operator W on the mode from j to I, W changes into 



W""^ = E c^„.i...M...«r'«r ■ ■ ■ a]_,''-'a,.r''^^^apa^^ ■ ■ ■ a^'af + h.c. (4) 
More concisely, it is 
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=: f{a{, ai, 4, 02 ■ ■ ■ a]_i, aj-i, aj, a] ■ ■ ■ ai, a}) : (5) 

The notation : ■ ■ ■ : indicates putting all the creation operators at left and annihilation 
operators at right for each mode in the functional /. 

Non-Gaussian states. The two mode squeezed state is a typical continuous variable state 
produced in laboratories.. Practically speaking, we can rarely have a pure squeezed state 
due to the noise. We formulate the following impure squeezed state from a noisy source : 

Mi,M2 

P= p{r) S (r) \00) {00\ S^r) dr + ^ p'{ni,n2)\ni,n2){ni,n2\. (6) 

ni,n2=0 

Here p{r) and p'{ni,n2) are the probabilistic distribution functions satisfying J p{r)dr + 
Snilnt^oP'(^i5 '^s) = 1, squeezing operator S{r) = exp[r(a|a2 — 0102)] {r > 0), tq and ri 
are positive numbers(ri > tq). Using the normally ordered form of squeezing operator we 
can simplify the integration part. Using Eq. we may perform the PT operation and we 
obtain 

. 1 Mi,M2 

pPT = p^r) ^ ^ tanh"i+"Var'ar|00)(00|afa^Mr+ ^ p'(ni, ^2)!^!, n2)(n2, ni|. (7) 

J cosh r ni,n2 ni,n2=0 

Consider the following linear superposed state in Fock space 

\^) = -L(|M,M + 2) - |M + 2,M)), (8) 
and M > Ml, M > M2. We have 

p{r 
'ro cosher 



(^Ip^^I^) = - r ^^te^nh^^'+'rdr < 0. (9) 

Jrn cosh r 



Thus state defined in Eq. is inseparable. As far as we have known, Eq. is a new 
formula for the impure squeezed state. This state is somewhat similar to the Werner state 
T7| for the discrete system. However, the result is different. A mixture of Bell state can 



be separable very often, while a mixture of different squeezed vacuum and number state is 
always inseparable. 

Besides the above noise due to the source, there are also noises due to the channel. We 
take the phase damping channel as an example here. Suppose initially we have a perfect 
pure squeezed state S'|00)(00|S'^ After solving the Master equation one obtains [|18| 



Pi't) = — 72— y tanh"i+"Vexp(-7t|rai ~ n2\'^)\ni,ni){n2,n2\. (10) 
cosh r„^^„2 

The entanglement formation of this p{t) has been numerically calculated However, as 
mentioned in [T^, "it is not clear from the present numerical analysis whether the state is 



3 



always entangled for finite 7^". The issue can be easily resolved in Fock space. According 
to Ref. |T^, if operator Q = tr/^p ® / — p is not positive definite, there is always a scheme to 



distill state p. We find here tr/)p(t) = ^^J^^ ^ J2n tanh "|n)(n|. It is easy to see (\l/2|f2|\E'2) < 0, 
taking |\E'2) = |(|00) + Thus we draw the conclusion: A two mode squeezed state is 

always distillahle( and inseparable) in phase damping channel. 

Gaussian states. We define the Gaussian state in the following normally ordered form in 
Fock space. 

"1 
2 



p =: exp 



{aRa + aBa^ + at^^a + ati?^c 



(11) 



a\,a\ - ■ ■ a\) 



Here a and a"^ are row vectors with / components, e.g. a = (ai, 02 ■ ■ ■ a^j; a- — ^u-i, ^2 
while a and are I component column vectors. Given a Gaussian operator that is not in the 
normally ordered form(e.g., the squeezed thermal state I^O])' always able to transform 

it into normally ordered form. There are a lot of well developed methods to do so(see e.g. 
|2T|,g). 

According to Eq. (|^), the transposition operation to the particles with subscript from j 
to / changes the state into 



=: exp 



1 ~ / R' B' 

- a, aT ^ 

2 ^ ' '\B'^ i?'t 




(12) 



where 





(13) 



and cr is a 2/ X 2/ matrix defined as u/f = 1, ai^f^i^f = 1 for / < j; (Jf^i+f = 1, cr/+/j = 1 
for / > j; and all the other matrix elements are 0. With the above formulae, one can take 
PT operation to any Gaussian state very easily. To check the inseparability, we need check 
the positivity of the Gaussian operator p^^ and we have 

Proposition : For the Gaussian operator p^^ defined above, it is positive definite if and 
only if all the eigenvalues of matrix B' not less than —1. 

Proof: If we find certain operator A so that p^^ = AA"^ , then p^^ is positive definite. Note 
that in the normally ordered form, operator : exp {aB'^a^ + aWa)/2 :=: exp[a'l"i?'a] :. 
Furthermore, if 5' > —1, we can use the formula 

: exp[(aS'^a^ + ^5'S)/2] := exp[^^t ln(5' + 

Denoting A = exp a^R'^ay2 exp ^a^ \ia{B' + l)a we have p^^ = AA^ . If some of the 
eigenvalues in matrix B' are equal to —1, considering the fact that : exp(— a|aj) := |0j)(0j|( 
|0j) is the vacuum state of the ith mode) and \/|0j)(0j| = |0j)(0j|, we can also easily factorize 
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into product form of certain A and A^. Suppose the first h eigenvalues of matrix B' are 
— 1, then we have 



A 



exp 



m'^ay2\ U : exp[(atA5)/2] :, 



(14) 



where An = —1 ii i < h, An > —1 if i > h, and unitary operator U satisfies U : exp[aAa''' + 
a^A'a] : W =: exp[ai?'^a^ + a'^B'a] :. The unitary operator U can be constructed in the 



following way: 

Suppose matrix B' can be diagonalized by unitary matrix D, i.e. B' = DAD^^. 
U to be the rotating operator satisfying Ua^W = a^D = 6t. Then we have 

: exp[^5'S] := E E T^^f^i = U ■ exp(^AS) : Ul 

k=0 i 



We define 



Using the fact y : exp(atAa :) =: exp(atAa/2) : we obtain Eq. (0). 

In the case that matrix B' has any eigenvalue less than —1, without loss of any generality, 
we suppose An < —1. Then the expectation value on state \ip) = exp[—aR'a/2]U\l, 0, ■ ■ ■ 0) 
is 



IP 



PTi 



;i| :exp[AnaIai] : |1) • (0| exp[A224a2]|0) ■ • ■ (0| exp[A,iafaz]|0). (15) 



Taking Taylor expansion of the exponential we know the terms of expectation values on 
vacuum state is 1, while the expectation value on state |1) is 



1 2 

;i| : exp[Anaiai] : |1) = (1|(1 + Aiia\ai + -Al^a\ al 



■)|1) = 1 + An <0. 



(16) 



So operator p^^ is not positive definite in this case. 

With the proposition one can easily check the inseparability of any Gaussian state. In 
general, one can easily obtain a sufficient condition for the inseparability of arbitrary state 
with them. In particular, it has been shown |]T^ recently that for N^l Gaussian state, there 
is no inseparable state that is positive definite after the PT operation. That is to say, our 
proposition have actually given the necessary and sufficient condition for the inseparability 
of bipartite Gaussian state which has N particles in one subspace and 1 particle in the 
other subspace. Moreover, connected with theorem 1 in Ref. |jl5|, the proposition gives the 
necessary and sufficient condition for the distillability for all bipartite Gaussian states. 
To demonstrate the calculation, we consider an un-normalized 3-variable Gaussian state 
defined as: 



p =: exp 



i<j 



Agajas 



(17) 



and 0<A<1,0<A3<1. If p is a physics state there are some restrictions to the parame- 
ters. Here we assume the parameters in p satisfy all those restrictions. Suppose mode 1 and 
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mode 2 are in subspace A while mode 3 are in subspace B. Taking PT operation to mode 3 
one easily obtains the matrix B' as defined by equation(0). The smallest eigenvalue (u;) of 
matrix B' for this state can be caculated easily. The condition for u < —1 is 



|Cl3|' + |C23r>(A + l)(A3+l). 



(18) 



Our proposition here is presented in Fock space. It can be easily applied to a type of more 
general problem, the entanglement property in a subspace of a global Gaussian state. The 
application background for this question is broad. For example, a multi particle bipartite 
state is prepared and then some of the particles are lost due to the operational errors, or 
some of the particles there are not accessible due to any other causes. We want to know 
whether the state in the accessible subspace is still distillable. For another example, the 
initially prepared Gaussian state interacts with its environment, a heat bath of harmonic 
oscillators. Finally we have a Gaussian state in the total system. But only a subspace is 
accessible. It is interesting to judge whether the state in the accessible subspace is distillable. 
The whole task is easy with our proposition in Fock space, because it is easy to obtain the 
explicit formula for a subspace state in a global Gaussian state, provided the global Gaussian 
state is expressed in the normally ordered form in Fock space. 

Bipartite state in a subspace of a global Gaussian state An /'—mode global 
Gaussian state pen can be defined in the form of equation (pT]). Modes numbered from 1 to 
/ belong to subspace C( which consists of subspaces A and B), modes numbered from / + 1 
to /' belong to the subspace D. We have 

1 



c 



PC = tiBpcD = tro : exp - [PcMcPc + PdMdPd + WoMof^c 
(ai, 02 ■ ■ ■ ctz, a|, 02 ■ ■ ■ aj), /3d = (a^+i, 0^+2 ■ ■ ■ «J+2 ' ' 



(19) 



''V )i 



C,D,0 Bc,D,0 
pt 

^C,D,0 ^C,D,0 



Mc,D,0 

, matrix Rc{Bc) is in / x / form, its matrix element Rcxx'{Bcxx') 



Rxx'{Bxx'),^ < x,x' < I, matrix R]j{Bu) is in (/' — I) x (/' — /) form, its matrix element 
Royy'iBDyy') = Ri+y,i+y' {Bi+y^i+y>) , I <y,y' < {I'-l), maXiix Ro{Bo) is in the form {l'-l)xl, 



its matrix element is Roz,z'{Boz 



Ro 



,l+z' 



),1 < z <l,l < z' <r -I. To take the 



partial trace in subspace D, we introduce the /' — / mode coherent state \Zd > in subspace 
D, \Zd >= \zi+i, zi+2 ■■ - zv >, ai+y\ZD >= zi+y\ZD >■ We have 



PC 



exp 



f3cMcf3c + Z*i,MDZn + 2(3cM^Zn\ Y[ d^zt/n : 

t=j+i 



(20) 



Using the complex Gaussian integration formula we obtain( up to a constant factor 
number): 



PC 



exp- 



f3c{Mc + 2M^M^'Mo)f3c 



(21) 
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We can easily recast equation (pT]) into the form of equation (0) through 




1 



) 



So we have given the exphcit Gaussian formula for the state in a subspace of a global 
Gaussian state. This equation, together with the known criterions on distillability and 
inseparability of Gaussian states gives the correspondings criterions in a subspace of a global 
Gaussian state. 

Concluding remark. In summary, we have presented a general method to detect the 
inseparability and distillability for states in Fock space. We have shown the inseparability 
for some non-Gaussian states. We give the necessary and sufficient condition for positivity 
of the normally ordered multi-mode Gaussian operator. This is also the necessary and 
sufficient condition for the distillability of all bipartite Gaussian states and the necessary 
and sufficient condition for the inseparability of (S> 1 Gaussian states. With the explicit 
formula for the bipartite state in a subspace of a Gaussian state given by us, all these 
criterions are applicable to a subspace inside the Gaussian state. 
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